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We extend the density-functional theory for superconductors (SCDFT) to take account of the 
dynamical structure of the screened Coulomb interaction. We construct an exchange-correlation 
kernel in the SCDFT gap equation on the basis of the random-phase approximation, where electronic 
collective excitations such as plasmons are properly treated. Through an application to fee lithium 
under high pressures, we demonstrate that our new kernel gives higher transition temperatures 
(T c ) when the plasmon and phonon cooperatively mediate pairing and it improves the agreement 
between the calculated and experimentally observed T c . The present formalism opens the door to 
non-empirical studies on unconventional electron mechanisms of superconductivity based on density 
functional theory. 

PACS numbers: 74.20.Mn, 74.25. Jb, 74.25.Kc, 74.62.Fj 



-Introduction. Electron mechanism of superconductiv- 
ity has long been a central topic of theoretical physics. 
Since the seminal discovery of cuprate superconductors, 
correlated systems with short-range Coulomb interaction 
have been extensively studied. On the other hand, sys- 
tems with long-range Coulomb interaction has also a long 
history, dating back to Kohn and Luttinger [T. Par- 
ticularly, there have been several theoretical proposals 
of superconducting mechanisms that exploit dynamical 
structure of the screened Coulomb interaction, which is 
represented by the frequency dependence of the dielectric 
function, e(u>). In fact, even without the electron-phonon 
interactions, it has been shown that superconductivity is 
induced by the dynamical Coulomb interactions due to 
exchange of plasmons [51 [3J or excitons [3]- More in- 
terestingly, the dynamical electron-electron interaction 
can cooperate with the ordinary electron-phonon pair- 
ing interaction. This possibility has been investigated 
for SrTiOa [5J [BJ, s-d transition metals [7J, insulator- 
sandwiched systems [H [5], and layered systems where 
soft electronic collective modes exist 10J. Moreover, re- 
cent discoveries of high-temperature superconductivity 
in doped band insulators have indeed stimulated further 
studies on this issue [TTtiM] . 

One important aspect of the above issues is that one 
must deal with quantitative competition between the 
attractive and repulsive contributions of the screened 
Coulomb interaction. Here, first-principles studies based 
on density functional theory are expected to provide 
deep insight. Recently, the progress in density functional 
theory for superconductors (SCDFT) [15l [16] has en- 
abled us to calculate superconducting transition temper- 
ature (T c ) nonempirically. Although the current SCDFT 
treats the screened electron-electron interaction within 
the static level, it achieves extreme accuracy for vari- 
ous conventional superconductors [TBTiLS] by including 
the dynamical phonon effect referring to the Migdal- 
Eliashberg (ME) theory [TS]. However, to formulate a 



quantitative ab initio method for unconventional plas- 
mon or exciton mechanism, we need to go beyond the 
currect SCDFT with the static approximation for the 
screened electron-electron interaction. 

In this paper, we extend the current SCDFT and de- 
velop a scheme to calculate T c with considering the dy- 
namical structure of the screened Coulomb interaction. 
Our formulation provides a general treatment of col- 
lective excitations described by the zeros of e(oj), that 
is, plasmons in solids. Our scheme enables stable cal- 
culations considering the frequency dependence of the 
screened Coulomb interaction, and is applicable to su- 
perconductors regardless of whether the plasmon effect is 
significant or not. As a representative system, we study 
superconducting fee lithium under high pressure [201 - 
[22] . which is expected to be similar to the electron 
gas with a small amount of carriers studied in Refs. [5] 
and [31 The previous calculation employing the current 
SCDFT [53] reproduced the experimental T c with the cal- 
culated electron-phonon coupling A>1 under the pressure 
>20GPa, whereas the most recent ab initio phonon cal- 
culation [53] gave far smaller A<0.8. We show that, even 
with A<0.8, the T c calculated with the plasmon contri- 
bution agrees well with the experimentally observed T c . 

-Formulation. In the current SCDFT QH HE], T c is 
obtained by solving the gap equation 

A 7 A 1 X-1C tanh[(/3/2)£W] 

n'VL' 

Here, n and k denote the band index and crystal mo- 
mentum, respectively, A is the gap function, and /3 is 
the inverse temperature. The energy E n k is defined as 
^«k=v / ^ k + A L and £«k = e«k - M is the one-electron 
energy measured from the chemical potential /i, where 
e„k is obtained by solving the normal Kohn-Sham equa- 
tion in density functional theory ^Ksl^nk) = Enkl^nk) 
with Hks and |y n k) being the Kohn-Sham Hamiltonian 
and the Bloch state, respectively. The functions Z and JC 
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FIG. 1: (Color online) Diagram corresponding to the electron 
nondiagonal kernel, JC cl . The solid line with arrows running 
in the opposite direction denotes the electronic anomalous 
propagator [15] . The red wavy line denotes the screened elec- 
tronic Coulomb interaction, which is a product of the inverse 
dielectric function e -1 and the bare Coulomb interaction V. 



are the exchange-correlation kernels describing the effects 
of the interactions. The nondiagonal kernel K. is com- 
posed of two parts /C=/C ph +/C el representing the electron- 
phonon and electron-electron interactions, whereas the 
diagonal kernel Z consists of one contribution 2=2^ 
representing the mass renormalization of the normal- 
state bandstructure due to the electron-phonon coupling. 
The phonon parts, K, ph and Z ph , have been formulated 
so that the conventional strong-coupling superconductiv- 
ity can be properly treated. The electronic nondiago- 
nal kernel JC cl corresponds to the diagram depicted in 
Fig. IT] with the static approximation e(oj)~e(0). We go 
beyond this approximation by retaining its frequency de- 
pendence. The diagram thus yields a new form for /C el 
as 



cl,dyn 
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where F„ k (iu;) = iu ^ Enk - iu> } Enk and w i and w 2 
denote the fermionic Matsubara frequency. Func- 
tion T4^„kr l 'k'(iw)=(^„ktV ;, n-ki|£ _1 (iw)V r |(p n 'k't < / 3 n'-k'i) 

denotes the screened Coulomb interaction scattering the 
Cooper pairs. 

In order to treat the plasmon-induced ui dependence 
of the screened Coulomb interaction, we employ the 
random-phase approximation (RPA [25 ) for e{iS). Fur- 
ther, we analytically carry out the Matsubara summa- 
tion in Eq. |2]) by approximating Wnkn'k' (iw) as a simple 
analytic function W^kn'k'Ci^)- We employ a multi-pole 
plasmon approximation represented as 



Wnkn'k' (hVm) = Wnkn'k' (0) 
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where N p is the number of plasmon modes and v m de- 
notes the bosonic Matsubara frequency. The functions 
&i;nkn'k' and Wi ; „kn'k' are the plasmon coupling and 
the plasmon frequency of the i-th mode, respectively. 
We have formulated so that W / n kn'k'(0)=Wnkn'k'(0)- 



Substituting Eq. ([3| in Eq. (pi), we finally obtain 
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where the function / is defined by Eq. (55) in Ref. [T3J 

We determine (aj ;n kra'k'i w i;nkn'k') in the following pro- 
cedure: (i) calculate the screened Coulomb interaction 
for the real frequency grid Wnk^'k' (w+i0 + ), (ii) determine 
the plasmon frequencies {wi; nk „' k <} by the position of 
the peaks up to the 7V p -th largest in ImW„kn'k' ( w +i0 + ), 
(ii) calculate the screened Coulomb interaction for the 
imaginary frequency grid VF„kn'k'(iw)) and (iv) using the 
calculated W n kn'W (iw), determine the plasmon coupling 
coefficients {ai ;7lk n'k'} via the least squares fitting by 
Wnkn'w (ico) . The approximate function Wnkn'k'(iw) de- 
rived in this procedure, as well as the analytically con- 
tinued ReIF„kn'k'( w + i0 + )> wen reproduces the calcu- 
lated W„kn'k'(iw) and ReWnkn'k'^ + iO" 1 "), respectively 
(see Fig. [2]). For numerical stability, we have imposed 
a constraint ai ;nkrl / k />0 and fixed ai ;nkrl < k ' to for the 
plasmon modes with peaks much lower than the high- 
est peak for each nkn'k'. We have considered up to two 
plasmon modes by setting N p =2 [26U29] , 
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FIG. 2: (Color online) Screened Coulomb interaction VF n kn'k' 
and the corresponding approximate function Wnkn'w f° r fee 
lithium in 14GPa calculated along the real frequency axis 
[(a)], and the imaginary frequency axis [(b)]. The band in- 
dices n and n' correspond to the partially occupied band, 
k=(2ir/a)(l/7, 1/7, 1/7), and k'=(0, 0, 0). 



-Application. We performed the calculations for fee 
lithium under pressure of 14, 20, 25, and 30GPa. All of 
our calculations were performed within the local-density 
approximation (301 131) using ab initio plane-wave pseu- 
dopotential calculation codes Quantum Espresso [521 - 
[33]. The lattice parameter used for each pressure 
was derived from the Murnaghan equation of state 
with total energy calculations for different lattice pa- 
rameters. Phonon frequency and electron-phonon cou- 
pling were calculated by the density functional pertur- 
bation theory |39j . The dielectric functions consider- 
ing the electronic screening were calculated within the 
RPA. The phonon contributions of the SCDFT exchange- 
correlation kernels (fC ph and Z ph ) were calculated us- 
ing the energy-averaged approximation |16j , whereas the 
electron contributions (/C cl ' stat and A/C cl ) were calculated 
by Eq. (13) in Ref.gOland Eq. @. The SCDFT gap equa- 
tion was solved with a random sampling scheme given 
in Ref. 1411 with which the sampling error in the calcu- 
lated T c was not more than a few percent. We also per- 
formed the calculation for aluminum in ambient pressure 
for comparison. 

Table |l] lists our calculated T c with £=/C ph (T c ph ), 
/C=/C ph +/C ol ' stat (T c stat ), and /C=/C ph +/C cl < stat +A/C cl 

N —2 

(T c p ). For examining the numerical stability of the 
fitting, the T c values calculated with A^ p =l (T c p_1 ) are 

iV —2 

also shown and they arc almost the same as the T c p 
values. We also show the electron-phonon coupling coef- 
ficient A, the logarithmic average of phonon frequencies 
win, the density parameter r s , and typical plasma fre- 
quency fip. Our calculated A well agrees with the latest 
calculation [2H, where the deviation is less than 0.04. 



TABLE I: Our calculated T c considering only the phonon 
contributions (ph) to the exchange-correlation kernels, 
the phonon and static electron contributions (stat), and 



all the contributions (N p 

,2 



2 and N p =l), Parameters 
exp[^ J dwlnuia 2 F(uj)/ui] were 



A=2 J cLua F(lj)/uj and win 

calculated from the Eliashberg functions a 2 F [15], and 
r s = -y/3/(47rp) and Qp — y^ATvp/m* were calculated using the 
electron density p and the band effective mass m* , with m* 
evaluated from the fitting of the calculated DOS by that of 
the parabolic band. 





Al 




fee Li 








14GPa 20GPa 25GPa 30GPa 


A 


0.417 


0.522 


0.623 


0.722 


0.812 


Win [K] 


314 


317 


316 


308 


304 




2.03 


2.71 


2.64 


2.59 


2.55 


fip[eV] 


16.2 


8.23 


8.44 


8.51 


8.58 


T ph [K] 


5.9 


10.0 


15.2 


19.0 


23.3 


ystat 

lf p=a [K] 


0.8 


0.7 


1.8 


3.2 


5.0 


1.4 


2.2 


4.4 


6.8 


9.1 


T c JVp=1 [K] 


1.4 


2.2 


4.1 


6.5 


9.1 


ycxpt. j K j 


1.20 a 


<4 




5-17 





15 







' C 

cxpts. 



□[20] 0(21] 
A[22] 0[43] <§ A 
A A 

o mo 



00 A 



13 



.Mi 



'Ref. 



20 25 
Pressure [GPaJ 

FIG. 3: (Color online) Our calculated T c (solid squares and 
circles) for fee lithium under high pressures compared with the 
experimentally observed values. The open symbols represent 
the experiments: Ref. [20] (open squares), Ref. [21] (open cir- 
cles), Ref. [55] (open triangles), and Ref. 43 (open diamonds). 



The material and pressure dependence of the experimen- 
tally observed T c is already reproduced in the calcula- 
tion with IC=JC ph , which can be understood in terms of 
the dependence of A. By including /C cl,stat , T c signifi- 
cantly decreases, and T c increases again when A/C cl is 
introduced. The origin of the systematic increase due to 
A/C el is discussed in the following. The marked difference 
of the dynamical contribution between Al and Li is sim- 
ply explained by the difference of f2 p : The characteristic 
structure in A/C el (given later) becomes relevant as Q p is 
lowered toward the phonon energy scale. 

We show in Fig. [3] the calculated T c for fee Li to- 
gether with the experimental values. By including A/C cl , 
the T c increases up to the range of experimental values 
(from solid square to solid circle). In contrast, in Al 
the agreement with the experiment is well achieved with 
/C=/C ph +/C cl ' stat , and AK, cl is irrelevant (see Table I). 
These results indicate the following: First, the high T c in 
fee Li under pressure is due to the help of the dynamical 
screening induced by plasmon, and second, our scheme 
works successfully regardless of whether their dynamical 
effects are strong or weak. 

In order to elucidate the origin of the increase of T c by 
the dynamical effect, we examined energy- averaged non- 
diagonal kernels Knk(g) = j^ J2n>W ^nkn'k"5(C - fn'k')> 

where N(£) is the electron density of states. With nk 
chosen as a certain point near the Fermi energy, we plot- 
ted the averaged kernel with A^ p =2 in Fig. El (a). We 
also decomposed the total kernel into /C ph , ^ sl . stat ) and 
A/C cl . Within the phonon energy scale, the total kernel 
becomes slightly negative due to /C ph , whereas it becomes 
positive out of the phonon energy scale mainly because of 
^el.stat The A/C el value is positive definite, but nearly 
zero for a low energy scale. As a result, a valley-like 
structure in the total kernel around the Fermi energy 
is enhanced. Thanks to the retardation effect [32] > this 
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FIG. 4: (Color online) (a) Contributions to the nondiagonal 
exchange-correlation kernel /C n kn'k' at T=0.01K calculated 
for fee lithium under pressure of 14GPa, averaged by equal- 
energy surfaces for n'k'. (b) The corresponding gap function 
calculated with (darker) and without (lighter) AfC cl . 



enhancement yields suppression of the effective repul- 
sion between quasiparticles. Remarkably, the low-energy 
valley-like structure in A/C cl appears in the energy range 
far smaller than the typical plasmon frequency (see Table 
[I]), which is consistent with the case of 3D electron gas 
reported by Takada [2]. 

The enhanced retardation effect depicted above can be 
seen more clearly from the gap functions plotted in Fig. [4] 
(b). In the SCDFT, the retardation effect is represented 
by sign inversion of the gap function in the energy region 
where /C„k(£) becomes positive [16]. The large negative 
gap value in the high energy region (darker points, >lcV) 
shows a significance of the retardation effect from this 
region. 

-Summary and conclusion. We extended the SCDFT 
to the dynamical-RPA level so that the plasmons in 
solids are considered. Stable T c calculations based on the 
present scheme were performed for fee Li, with which the 
plasmon effect on T c was quantified. We found that this 
effect substantially raises T c by cooperating with phonon. 
In addition, we showed that the calculated T c agrees bet- 



ter with the experimental values than that within the 
static-RPA level. Here, the plasmon-induced dynamical 
component of the screened Coulomb interaction enhances 
the retardation effect. Our present scheme provides a 
firm foundation for density functional theory for uncon- 
ventional superconductivity induced/assisted by the elec- 
tron correlation. 
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